Renormalization of a gapless Hartree-Fock approximation to a theory with 

spontaneously broken 0(7V)-symmetry 
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The renormalization of a gapless "I>-derivable Hartree-Fock approximation to the 0(iV)-symmetric 
\<f) 4 theory is considered in the spontaneously broken phase. This kind of approach was proposed by 
three of us in a previous paper Q| in order to preserve all the desirable features of "I>-derivable Dyson- 
Schwinger resummation schemes (i.e., validity of conservation laws and thermodynamic consistency) 
while simultaneously restoring the Nambu-Goldstone theorem in the broken phase. It is shown that 
unlike for the conventional Hartree-Fock approximation this approach allows for a scale-independent 
renormalization in the vacuum. However, the scale dependence still persists at finite temperatures. 
Various branches of the solution are studied. The occurrence of a limiting temperature inherent in 
the renormalized Hartree-Fock approximation at fixed renormalization scale [i is discussed. 
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I. INTRODUCTION 

Self-consistent "^-derivable approximations were intro- 
duced long ago in the context of the nonrelativistic many- 
body problem 0, Q and then extended to relativistic 
quantum field theory 0, . Recently the interest in this 
method has been revived in view of its fruitful appli- 
cations to calculations of the thermodynamic properties 
of the quark-gluon plasma lol and to non-equilibrium 
quantum- field dynamics 0, H- El ■ i n particular in terms 
of the off-shell kinetic equation |lCj . 

^-derivable approximations are preferable for the dy- 
namical treatment of a system, since they fulfill the con- 
servation laws of energy, momentum, and charge 0, 0, 
ITol ]. Moreover, the <f>-derivable scheme also guarantees 
the thermodynamic consistency of an approximation [j| , 
which makes it advantageous also for thermodynamic 
calculations. However, the <I>-derivable scheme has its 
generic problems which were also realized long ago pSHTlj. 

The first problem is related to the fact that ^-derivable 
Dyson-Schwingcr resummation schemes violate Ward- 
Takahashi identities beyond the one-point level. This, 
in particular, results in the violation of the Nambu 
Goldstone (NG) theorem % GH E in the phase 
of spontaneously broken symmetry. On the other hand, 
so called "gapless" approximations respect the NG 
theorem (which is referred to as the Hugenholtz-Pines 
theorem in physics of Bose-Einstein condensed systems), 
though violate conservation laws and thermodynamic 
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consistency. In Ref. 0| it was shown that any 
derivable approximation can be corrected in such a way 
that it respects the NG theorem and becomes gapless. 
However, such modifications again violate conservation 
laws and thermodynamic consistency and, hence, leads 
back to the problems of the gapless scheme. Recently 
three of us have proposed a phenomenological way to 
construct a "gapless <I>-derivable" Hartree-Fock (gHF) 
approximation to the A0 4 theory in the phase of spon- 
taneously broken O(N) symmetry This approxima- 
tion simultaneously preserves all the desirable features of 
<I>-derivable schemes and respects the Nambu-Goldstone 
theorem in the broken phase. The treatment of Ref. pj 
was based on a naive renormalization, where all diver- 
gent terms were simply omitted. This was done in order 
to avoid possible confusions between effects of restoring 
the NG theorem and those related to renormalization. In 
the present paper we return to the issue of renormaliza- 
tion. 

The renormalization of the ^-derivable approximations 
is precisely the second main problem. Following Baym 
and Grinstein || it was believed that renormalization of 
<I>-derivable approximations is possible only with medium 
(e.g., temperature) dependent counter terms, which is 
inconsistent with the goal of renormalization. Great 
progress in the proper renormalization of such schemes 
was recently achieved in Refs. [II EH ES ■ As 
the main result it was shown that partial resumma- 
tion schemes can indeed be renormalized with medium- 
independent counter terms provided the scheme is gener- 
ated from a two-particle irreducible (2PI) functional, i.e., 
a ^-functional. Still, as we are going to demonstrate be- 
low, certain problems remain in the case of spontaneously 
broken symmetry. 

As an example case we investigate the O(N) model 
in the spontaneously broken phase which is a traditional 
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touchstone for new theoretical approaches, well applied 
to a variety of physical phenomena, such as the chiral 
phase transition in nuclear matter. Thus we continue 
the discussion of the "gapless Hartree-Fock approxima- 
tion started in the previous paper Q] and investigate its 
features towards renormalization in comparison to the 
standard HF-approximation. 



II. GAPLESS HARTREE-FOCK (gHF) 
APPROXIMATION 



We consider the 0(N)-mode\ Lagrangian 
1 
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where <fi = 

A2 _ ± ± 



>i,4>2, 4>n) is an TV-component scalar field, 
with summation over a implied. For H = 
this Lagrangian is invariant under O(N) rotations of 
the fields. If H = and m 2 < 0, the symmetry of the 
ground state is spontaneously broken down to 0(N — 1), 
with N — 1 Goldstone bosons (pions). The external field 
H ■ 4> = H a <fi a is a term which explicitly breaks the O(N) 
symmetry. It is introduced to give the physical value of 
140 MeV to the pion mass. 

The effective action T for this Lagrangian in the real- 
time formalism is defined as (cf. Ref. Q) 



r^Gj^/oW + ^TrtlnG- 1 ) 

+ ^Tr (D~ l G - 1) + $ rea l-time{^ Cf), 



(2) 



where 4> is the expectation value of the field, G is the 
Green's function, D is the free Green's function, Io(4>) 
is the free classical action of the <fi field, Tr implies 
space-time integration and summation over field indices 
a,b, ... All the considerations below are performed in 
terms of the thermodynamic $ functional which differs 
from ^cai-timc in the factor of i/3, where (3 = 1/T is 
the inverse temperature. In the case of a spatially ho- 
mogeneous thermodynamic system, an additional factor 
appears: the volume V of the system. Thus, the thermo- 
dynamic $ is 



$ = (-iT/F)$ real - 

In terms of the CJT formalism [3 
action T is given as (e.g., cf. Ref. Q) 

r{0, G} = /(<£) + ^Tr (In G- 1 ; 



the same effective 



^(d-'g-i) +r 2 {0,G}, 



(3) 



i.e., already in terms of the tree-level Green's function, 
Dj,(l,2) = 5 2 I/54>20~(t>i, and the full classical action of 



the <j) field, I(<f>). In the thermodynamic limit, the ef- 
fective potential, Vqjt, and its interaction part, V 2 , are 
defined as 

^cjt = (-iT/V)T, V 2 = {-\T/V)T 2 . 

Naturally, V2 is similar to 3? but is not quite the same. 
Contrary to V 2 , the <3? functional includes all 2PI inter- 
action terms, i.e. also those of zero and first loop or- 
der which result from interactions with the classical field 
(first two graphs in 

The gHF approximation to the O(N) theory is defined 
by the <& functional Q 



$gHF = 



+ A$, 



(4) 



where the diagrams on the r.h.s. constitutes the con- 
ventional HF approximation, while the phenomenological 
NG -theorem-restoring correction A$ is specified below, 
see Eq. (JSJ. Here the crossed pins denote the classical 
fields (p a , and loops are tadpoles 
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d 4 kG ab (k) 



(5) 



in terms of G a b Green's functions, where the Matsubara 
summation 

f d\f(q)=T jr [^f(2rinT,q) (6) 

is implied with T being a finite temperature. 

Within the ^-derivable scheme the r.h.s. of the equa- 
tions of motion for the classical field (J) and the Green's 
function (self-energy S) follow from the functional vari- 
ation of $ g HF with respect to the classical field cj> and 
Green's function G, respectively 



D(j> + m 2 



G^-D- 1 = 




(7) 



(8) 



where D is the free propagator. 

The Ai> correction, introduced in Ref. Q, is unam- 
biguously determined proceeding from the following re- 
quirements: (i) it restores the NG theorem in the broken- 
symmetry phase, (ii) it does not change results in the 
phase of restored O(N) symmetry, because there is no 
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need for it, (iii) it docs not change the HF equation for 
the classical field, since the conventional ^-derivable and 
gapless schemes 0, El provide the same classical-field 
equation already without any modifications. In partic- 
ular, due to this latter requirement the A$ correction 
does not contribute to the classical-field equation 0. 
Proceeding from these requirements, this A$ can be pre- 
sented in manifestly O(N) symmetric form 

^ = -^[N(Q ab ) 2 -(Q aa ) 2 ]. (9) 

Here and below, summation over repeated indices 
a, 6, c, ... is implied, if it is not pointed out otherwise. 

The nature of this correction can be understood as fol- 
lows. For the full theory, i.e., when all diagrams in the 



$ functional are taken into account, the gapless and <&- 
derivable schemes are identical and both respect the NG 
theorem. The conventional ^-derivable HF approxima- 
tion omits an infinite set of diagrams which is necessary 
to restore its equivalence with the gapless scheme. The 
A4> correction to the HF approximation takes into ac- 
count a part of those omitted diagrams (at the level of 
the actual approximation), and thus restores this equiv- 
alence in the pion sector. For the further discussion we 
switch to the notation in terms of the CJT effective po- 
tential, see e.g. 0,0, in order to comply with pre- 
vious considerations in the literature. 

The manifestly symmetric form of the CJT effective 
potential in the gHF approximation reads 
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k2 f-H -<j)+- { ^fclndetG-^fc) 
2 J 8 

G ba (k) - 1 



(k 2 + m 2 ) 5 ab + — {<f 2 5 ah + 2(t) a (j) b ) 



+ (QaaQbb + ZQabQba) + A$, 



(10) 



e.g., cf. |20j| . All quantities are symmetric with respect 
to permutations of indices. 

Since the self-energy © is momentum independent, 
the general form of the Green's function can be written 
as follows 



G- b \k)=k 2 S ab + M 2 ab , 



(11) 



where M 2 h is a constant mass matrix. The equations 
for G c d; i.e., for the corresponding tadpoles Q ab , and the 
fields 4> c result from variations of V^hf over G cc i and <f> c , 
respectively, 

M 2 d = m 2 S cd (12) 

+ ^ W 6 cd + tycfa + SQaaScd + 2(1 - N)Q cd ] , 

H c = m 2 (j) c + A [0 2 c + Q aa (j) c + 2Q cd <jy d ] . (13) 

These are equations in a general nondiagonal representa- 
tion. Applying projectors 

Kd = j^j (Scd - Md/tf) , (14) 

n° d = c d /0 2 , (15) 

to Eq. l|T2f> . we project it on ir and a states. In order to 
project the mean-field equation (|13f) on the er-direction, 
we just multiply it by <j) c . 

In the diagonal representation (<j> a ^ 0, H a = H and 
H-n = 4> n = 0) these equations take the following form 

Ml = m 2 + A [30 2 + (5 - 2N)Q a + 3(N - l)Q ff ] 



M 2 + - [2<j> 2 + 2(N — l)(Q 7r — Qa)] , (16) 



M 2 = 
H = 



- [4> 2 + SQa + (N- 1)Q W ] 



N 



+ 3Q CT + (N- 1)QJ 



(17) 
(18) 



where M 2 = W dc Ml d and il/ 2 = W dc Ml d . Here we used 
Qa = Qaa and = Q VTC in terms of definition JSJ. 
From these equations it is evident that the NG theo- 
rem is fulfilled. Indeed, in the phase of spontaneously 
broken symmetry (H = 0) the square-bracketed term 
of the field equation (|18J) equals zero, which is precisely 
the pion mass, cf. Eq. (|17J) . At the same time, as it 
has been demonstrated in numerous papers (see, e.g., 
Refs. Hlll|ll|l3), the solution of the conventional HF 
set of equations 0-©, i.e., without A$, violates the 
NG theorem. A detailed analysis of the conventional HF 
equations in a notation similar to ours has been given in 
Ref. [H. 



III. RENORMALIZATION OF THE gHF 
APPROXIMATION 



Significant progress in proper renormalization of 
^-derivable approximations was recently achieved in 
Refs. 0, fill Il6j. Here we follow the renormalization 
scheme of Ref. [l3|], i.e., that constructed precisely for 
the conventional HF approximation to A(/> 4 theory in the 
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O(N) broken phase. This renormalization is based on 
the BPHZ formalism. 



A. Equations of Motion 

The equations of motion i|16|) - (|18[l involve tadpole 
terms which, based on the explicit form of the Green's 
function can be written as 



Qa 



d 3 k 



(2t) 3 e a (k) 



n(e a (k) 



where e a (k) = {P + Ml) 1 / 2 and 

1 



n(e) 



exp(e/T) - 1 



(19) 



(20) 



is the thermal occupation number. Evidently, the Q- 
function consists of two parts 



where 



Ql 



Qa = Qt+< 



d 3 k 



W e a (k) 



)(div) 



n (e a (k) 



(21) 



(22) 



is the convergent thermal part of the tadpole, which is 
finite, and the divergent part 



(dr 



i r d 3 k i 
2 J WfTJk) 

M 2 



(4tt) 2 



1 , Ml 



(23) 



is regularized within dimensional regularization. Here 
e — > 0, and /i is a regularization scale. We apply the 
same mass independent renormalization conditions as in 
Ref. |l3j. i.e., that in the symmetric vacuum the self- 
energies vanish [2l| 

T, a (T = 0,0 = 0,?ti 2 = fi 2 > 0) = 0, (24) 
d m 2E a (T = 0,0 = 0,m 2 = ii 2 > 0) = 0, (25) 

for all a. Such a renormalization scheme preserves the 
0(2V) symmetry of the model [2l]]. Since S a is momen- 
tum independent in the approximation under considera- 
tion, additional momentum-derivative conditions are not 
required. Upon application of this scheme, Eqs. I|16[) - 
(|T%|l keep their form with the Q a quantities substituted 
by the renormalizcd tadpoles. 



(rcn) 



(4tt) 2 



Mi In 



Ml 



-1 



(26) 



As it was shown in Ref. [l3[ . this renormalization de- 
scription requires only vacuum (temperature indepen- 
dent) counter terms. For the sake of further discussion, 
note that according to Ref. 01 the renormalization of the 
conventional HF approximation results in precisely the 
same equations as those in the CT scheme of Ref. [ljj , in 
spite of the different approaches used. The renormalizcd 
conventional HF approximation was thoroughly studied 
in |l9j |. Therefore, those results are very useful for com- 
parison with the present treatment. 



B. Effective Potential 

The thermodynamic potential l|10[) is renormalizcd fol- 
lowing the procedure outlined in Ref. 0] . In the gHF ap- 
proximation complications arise only in the lndet G~ 1 (k) 
term. Because of the topology of the diagrams used in 
the gHF approximation, for all other contributions to the 
effective action, wc only have to insert the already renor- 
malized self-energies, i.e., Q [ a cn) tadpoles of Eq. ftfy. 
to renormalize them. This is legitimate, because we do 
not encounter any additional contributions from subdi- 
vergences, cf. Ref. I n the diagonal representation 

the remaining part to be renormalized takes the form 

- / rf 4 fclndctG _1 (fc) = L a + (N - l)L n , (27) 
where we have introduced the brief notation 

= /(0{f+ ri °H*>>(-|)]}' < 28 > 

The L a also consists of two parts: the convergent thermal 
part 



(2tt) 3 



1 — exp 



e a (fc) 
T 



(29) 



which is finite, and the divergent integral 

r(div) _ f d 3 k e a (fc) 
(2tt) 3 2 ' 



(30) 



This expression implicitly depends on temperature 
through the mass M a . We regularize it by means of a 
momentum cut-off A: 



5 



8 



d 3 k e a {k) 
(2tt) 3 2 

2AVA 2 - M 2 - 5AMVA 2 - AP - 3M 4 In M + 3 A/ 4 In ( A + ^A 2 - A/ 2 



e(A 2 - £ 2 (A?))=^((A 2 -M 2 ) 3/2 A 



(31) 



To renormalize the effective potential, we use a mass-independent renormalization scheme in order to avoid effects of 
unphysical IR singularities. We impose the following renormalization conditions on the effective potential 



y (rc„) ( T = O ,0 = O,m 2 ) 
d m 2V {lcn) (T = O,0 = O,m 2 ) 
a 2 , 2 l/ (rcn) (T = O,0 = O,m 2 ) 



m 2 =^ 2 >0 


= o, 


m 2 =[i 2 >0 


= o, 


m 2 =n 2 >0 


= 0. 



(32) 
(33) 
(34) 



We need precisely these three conditions to remove all the divergences. Imposing these conditions, we keep in mind 
that other parts of the effective potential, except for L a , have already been renormalizcd such that they fulfill these 
conditions on their own. This leads to 



^(ren) 



L, 



lim 

A^oo 
1 

1287T 2 



4 res) (M Q ) - 4 rcg V) - {mi - ^ 2 )^ M - V M « ~ t ^wr* 

3M 4 - 4M V + M 4 - 2TV/ 4 In ^ 



(35) 



In terms of these La Cn ^ and Qa the renormalized effective potential reads 



T) = \m*cp + ^ -Hj> + l£»> + (TV - 1)4^) 

(Q£' cn) + (TV — l)Qt n) ) - MlQ^ — (TV — l)A/ 2 Qi rcn) + ^cf> 2 (sQ^ + (TV - 1)Q^ 
(Q(;' cn) + (TV — l)Qi rcn) ) 2 - 2(TV - 1) ([Q£ en) ] 2 + (TV - 1) 
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2 

A 
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)(ren) 



(36) 



C. Vacuum (T = 0) 



At T = 0, the quantities under investigation are "experimentally" known 1 : M n {T = 0) = m„ = 139 MeV, 
M a (T = 0) = m a = 600 MeV, and the pion decay constant <fio = f n = 93 MeV. At T = these known quantities 
should satisfy Eqs. (|16|H(|18[) with renormalized tadpoles Qa Bn \ cf. Eq. 126|) . 



H 

U 



2A 2 2(TV-1)A 
TV /7r (4tt) 2 TV 



In- 



A j 2 -| — ^ — 

TV (4tt) 2 TV 

Aj 2 -| — - — 

TV (4tt) 2 TV 



(/V + 2) M 2 + 3m 2 In 



(N + 2) M 2 + 3m 2 In 



m 

,,2 



/'- 



In order to make this set of equations consistent, we should put 

H = mlU, 



In- 



+ (TV - l)m 2 In 



+ (TV - l)m 2 In 



(37) 
(38) 
(39) 

(40) 



i.e., precisely the same as at the tree level. Then Eqs. (13811 and l|39|l become identical. Now we can solve for the 
remaining equations, (|37|l and l|38|) . and express the renormalized quantities m and A in terms of physical quantities 



These values arc relevant for the case N = 4. 
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m w , m<j, and f„ 



N (ml - ml) 2ft 



2{N - 1) 



A 



N (4tt) 2 7V 



(An) 2 
(N + 2) A i 2 



3m a hn^ _ 1 + (JV- l)ml ln^ - 1 



(41) 
(42) 



The expressions manifestly show that the performed 
renormalization is /^-scale independent in the vacuum. 
Indeed, m 2 and A can take any values depending on the 
choice of the renormalization scale fi, while the same ob- 
scrvables keep their physical values at any renormaliza- 
tion scale fj, (within a certain range of fi). In this respect, 
the gHF approximation is similar to the leading-order 
1/iV-approximation, where also both the NG theorem 
and scale-independent renormalization in the vacuum are 
fulfilled 0| . The restriction to a certain range is related 
to the conditions A > and m 2 < which should be met. 
At the scale /zq, determined by the equation 



J TV 



(N-l) 
(4tt) 2 



m„ In 



Mo 



2 1 ""<T 

— m In — s- 

Mo 



= 0, (43) 



A and m 2 become singular, cf. Eqs. (|41(l and (|42J) . More- 
over, /i < /io implies m 2 > and A < 0, which defers a 
spontaneously broken phase and makes the theory un- 
stable because of uncompensated attraction (A</> 4 < 0). 
Therefore, the range of scale independent renormaliza- 
tion in the vacuum is restricted from below by 



M > Mo- 



(44) 



The situation is completely different within the conven- 
tional HF approximation. In that case, physical obscrv- 
ables do depend on the renormalization scale and can be 
reproduced only at a single value of the scale fi 2 = ml/e, 
as demonstrated in [J^. The reason of this difference is 
that the conventional HF approximation of |19| violates 
the NG theorem at H = 0. This stems from the fact that 
the equations for the pion self-energy and the mean field 
are not identical in the broken phase and thus leave three 
nondegenerate equations for three quantities m 2 , A and 
(i, from which the \i value is unambiguously determined. 
In our case, this set of equations is degenerate, i.e., the 
equations for the pion self-energy and the mean field are 
identical as a consequence of the fulfilled NG theorem. 
This gives us freedom for an arbitrary choice of (i. Thus, 
we arrive at an important conclusion which concerns all 
partial resummation schemes applied to the case of spon- 
taneously broken symmetry: a scale-independent renor- 
malization in the vacuum is possible only if the scheme 
preserves the NG theorem. This is an important aspect 
with respect to possible renormalization-group consider- 
ations for this type of self consistent approximations. 



However, the scale dependence still persists at finite 
temperature, which is already seen from the analysis of 
the symmetry-restoration points. 



D. Symmetry Restoration Points at H — 

Starting from the broken phase, where M 2 = and 
<f) 2 > 0, there exists a temperature range T R G [T\, T 2 ], cf. 
Fig. 2] below, where the classical field vanishes together 
with the pion mass 



M 2 (T R ) = 4> 2 {T R ) = 0. 



(45) 



This precisely occurs, when the two equations l|17|) and 
(|T5|l reduce to a single one. Solving for this single equa- 
tion with renormalized Qa™^ tadpoles l|26|l gives 



T 2 R 



12 



(N + 2) 



3m| 
(4tt) 2 



In 
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3 M 2 (T R ) 
(N - 1) ml 
M 2 (T R ) 



J n 



(46) 



where we have used Eqs. I|41ll and l|42() to express m 2 and 
A in terms of the vacuum mass m a . Strictly speaking, 
this is not a solution for T R , since the r.h.s. of l|46|) 
still depends on T R through M 2 (T R ). Nevertheless, this 
expression is already quite simple to analyze. 

The lower bound T\ of condition (|45J) is of central im- 
portance for thephase transition in the conventional HF- 
approximation [l| but of minor relevance in the gHF- 
scheme, since it corresponds to the metastable solution, 
cf. Sect. |^ below. It is determined if simultaneously 
also M 2 = occurs 



12 



(7V + 2) 



fl 



3ml 
(4tt) 2 



In- 



(47) 



It is still /i-dependcnt, in spite of the scale- independent 
renormalization in the vacuum. At large fj,, i.e., above 
some fix, T 2 can even become negative, which means that 
then this solution does not exist. 

For the stable solution of the gHF approximation it 
numerically occurs that M a (T R ) pa m a , cf. Ref. pj and 
Sect. ITVl below. This almost removes the M,j(T R ) de- 
pendence from the r.h.s. of Eq. (|46|l and makes the 
corresponding temperature Ti almost /i-indepcndcnt 



To 



12 



(N-l) 



f 



(48) 
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This value coincides with that of the naive renormaliza- 
tion pj . The solution T 2 corresponds to a partial symme- 
try restoration, since here we still have M a (T 2 ) ^ M T (T 2 ) 
in spite of 4> = 0. 

IV. RESULTS FOR N = 4 

For the numerical calculations we use the following pa- 
rameters: m a = 600 MeV and f n = 93 MeV. The pion 
mass is either zero, m„ = 0, in the case of exact symme- 
try, or m T = 139 MeV for the approximate symmetry. 
The general structure of the solutions to the renormal- 
ized Eqs. H16f) - I|18[l is similar to that obtained with the 
naive renormalization Q but with extra complications 
caused by the additional dependence on the renormaliza- 
tion scale fj,. According to Eq. l|4^|) physically reasonable 
solutions exits only for /i > (Iq, where no as the solution 
of Eq. iflTSJl equals ~ 200 MeV for the above specified 
parameters. 

There are several different branches of the solution. 
Stable and physically meaningful are determined by the 
principle of maximum pressure, the pressure being given 
by the effective potential ijHB^I 

P = -V g ( Hp ) (0 1 T) + const. (49) 

Here the constant is determined by the condition that 
the pressure should vanish for the physical vacuum, i.e., 
in the spontaneously broken phase, while our renormal- 
ization condition l|32|) determines Vg^f* to be zero in the 
unphysical, symmetric vacuum. 

A. Exact 0(4) Symmetry 
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FIG. 1: Meson masses as functions of temperature for /i = 
600 MeV and rrw = case. Stable branch is presented by 
solid lines, whereas the metastable one - by dashed lines. 
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T [GeV] 

FIG. 2: The same as in Fig. Qbut for the pressure difference 
between stable and metastable branches. 



The actual structure of the solution depends on the 
renormalization scale \x. We start with a moderate scale 
fx = 600 MeV, i.e., of the order of m a . The results are 
presented in Figs. ^El 

In the narrow temperature range, displayed in Fig. ^ 
the results are qualitatively similar to those obtained 
with the naive renormalization Q]. The stable branch 
starts at T = from the physical vacuum values for the 
masses and the classical field and crosses the metastable 
branch at T cross ps 440 MeV. In terms of the pres- 
sure, they are touching rather than crossing (see Fig.|5J). 
Therefore, no transition from one branch to another oc- 
curs at T cross . In the broken-symmetry phase, the pion 
mass equals zero. Then a phase transition of the second 
order occurs at T 2 ~ 180 MeV, at which the field be- 
comes zero (see Fig. |3J). However, the ir and a masses 
still differ beyond this transition point. They become 
equal only after a second phase transition, which is also 
of second order, at T cross . Note that the equal-mass solu- 
tion above T cross is precisely the same as in the conven- 
tional HF approximation (cf. Rcfs. since the 
gapless modification term © vanishes in this case. The 



T\ point proves to be irrelevant for the stable branch. 
Rather it is the starting point for the metastable branch, 
which in the range of T± < T < T cross precisely coincides 
with the solution of the conventional HF approximation 
[l3l ITflj . The corresponding field is always zero for this 
branch. Contrary to the case of the naive renormaliza- 
tion Q , this metastable branch ends at some temperature 
(ps 650 MeV). 

However, in the wider temperature range, displayed in 
Fig. 01 we see a significant difference from the results 
obtained with the naive renormalization [JJ . Neither sta- 
ble nor metastable solutions exist above a certain tem- 
perature Tend which is w 1.15 GeV for the considered 
fi = 600 MeV. Let us remind that this solution above 
Tcross corresponds to the conventional HF approximation. 
Therefore, it is precisely the same as in previous conven- 
tional HF calculations with renormalization [HEil. The 
occurrence of such an end point in the HF approximation 
was first pointed out by Baym and Grinstein 

At T cn d, the stable branch of the solution joins the up- 
per branch. This upper branch at any temperature corre- 
sponds to equal masses and zero field, i.e., = M„ and 
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FIG. 3: The same as in Fig. [T]but for the field <f>. 
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FIG. 4: Meson masses of Fig. but in wider temperature 
region. The upper metastable branch is displayed by the long- 
dashed line. 
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FIG. 6: Meson masses as functions of temperature for /i = 
1200 MeV and m-^ = case. Stable branches are presented 
by solid lines, whereas metastable ones by dashed lines. 



= 0, and hence is also a solution to the conventional 
renormalized HF approximation. It starts with very high 
values of masses (« 3 GeV) in the vacuum and also ends 
at T cn( j. The vacuum pressure for this branch, evaluated 
according to Eq. (|4*Tj|) . is rather high (see Fig. [5j|. 

It is important note that for the upper branch and at 
temperatures near the endpoint the logarithmic terms in 
the gap equation, oc \n(M^ a /fi 2 ), become large. This in- 
dicates that at such points the expansion of the $ func- 
tional in powers of the renormalized coupling becomes 
unreliable, because the effective coupling becomes large. 
Therefore we consider the upper branch not a physically 
meaningful solution. The same holds true for tempera- 
tures close to the endpoint temperature. Such a behavior 
must be expected for any effective theory and was in- 
deed also observed in Quantum Hadro Dynamics (QHD) 
in 

At larger renormalization scales, [i, the global pattern 
of the solution remains qualitatively similar, as seen from 
Fig. [51 Only T cross and T cn( j move to higher temper- 
atures. Inspecting the low-temperature region in more 
detail, cf. Fig. [7| we see that a new metastable solution, 
which ends already at rather low temperature, appears. 



This new solution has a nonzero field (Fig. [SJ) but violates 
the NG theorem. In addition the metastable solution, 
which before started at T\, now begins at zero tempera- 
ture, since at /i = 1.2 GeV we have already T\ 2 < 0, cf. 
Eq. (gZJ. 

If we take a low value for the scale /i, the structure of 
the stable solution becomes more involved, see Figs.l^llTTl 
In the broken-symmetry sector we still have a massless 
pion and a nonzero field, see Fig. 1101 However, at higher 
temperatures the metastable branch, displayed by the 
dotted line, reveals back-bending. This back-bended part 
of the branch turns out to be the most stable one, see 
Fig. 1111 As a result we arrive at a complicated structure 
for the stable solution, where even the pressure turns out 
to be discontinuous. 

A common feature for all scales is that the point of the 
first phase transition, T2 ~ 180 MeV, is approximately 
/i-independent and has about the same value as that in 
the naive renormalization scheme At the same time, 
the point of the second phase transition, T cross , and the 
end point of the stable solution, T on( j, arc essentially /1- 
dependent. 
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FIG. 7: Zoomed low-temperature region of Fig. HU 
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B. Approximate 0(4) Symmetry 

In the case of explicitly broken symmetry (m^ = 
139 MeV), the structure of solutions at various /i scales is 
similar to that described above for the chiral limit. Even 
the behavior of metastable branches remains similar. We 
illustrate the changes on the example of fj, = 600 MeV, 
see Figs. I12H14I which looks the most physically appeal- 
ing and is close to results of the naive renormalization 
The main difference from the m T = case is that 
the sequence of two phase transitions is transformed here 
into a smooth cross-over transition. 



V. CONCLUSION 

We have studied a renormalized version of the gapless 
<I>-derivable HF approximation to the A</> 4 theory with 
spontaneous breaking of the O(N) symmetry, proposed 
in Ref . Q . This gHF approximation simultaneously pre- 
serves all the desirable features of <I>-derivable approx- 
imation schemes (i.e., the validity of conservation laws 
and thermodynamic consistency) and respects the NG 
theorem in the phase of spontaneously broken symmetry. 



FIG. 9: Meson masses as functions of temperature for /i = 
300 MeV and = case. Stable branches are presented by 
solid lines, whereas metastable ones - by dashed lines. 
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FIG. 10: The same as in Fig. ©but for the field 



This is achieved by adding a correction A<E> to the con- 
ventional $ functional. The nature of this correction can 
be understood as follows. The conventional ^-derivable 
HF approximation cuts off an infinite series of diagrams, 
among which are those providing the NG theorem in the 
phase of spontaneously broken symmetry. By introduc- 
ing the A<!> correction to the HF approximation we take 
into account a part of those omitted diagrams (at the 
level of the actual approximation), which restores the 
NG theorem. 

An advantage of the gHF approximation is that it al- 
lows for scale-independent renormalization in the vac- 
uum, unlike the conventional HF approximation |l9j| . 
The scale independence in the vacuum is a direct con- 
sequence of the NG theorem, which makes the equations 
for the pion self-energy and classical field degenerated. 
However, even in the gHF approximation, only renor- 
malization scales higher than a certain value (no, cf. Eq. 
l|4H|l) are allowed, in order to ensure stability of the renor- 
malized approximation. 

Nevertheless, the scale dependence still persists at fi- 
nite temperatures. The violation of renormalization-scale 
independence of <f>-dcrivablc approximations was shown 
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in \2?t l24j from the point of view of the renormalization- 
group [3 function. There the (3 function, evaluated from 
the ^-functional formalism, was shown to deviate from 
its perturbativc expansion, beginning at orders in the ex- 
pansion parameter, higher than that explicitly taken in 
the <E> functional. The reason is the violation of "crossing 
symmetry" in the sense of 01 : Solving the self-consistent 
equations of motion corresponds to a partial resumma- 
tion of diagrams to any order in the expansion parameter 
(e.g., the coupling constant A or ft, i.e., the order of loops 
in perturbative Feynman diagrams) which is necessarily 
incomplete for any truncation of the <J> functional. 

Within our renormalization scheme, it becomes clear 
that the renormalization-scalc dependence at finite tem- 
peratures originates from the subtraction of the "hidden 
subdivergence" of the four-point function inside the self- 
consistent tadpole loop. As shown in |14j . this four-point 
function consists of a resummation in only one channel, 
and thus the (3 function of this resummed four-point func- 
tion deviates from the correct one at orders higher than 
contained in the approximation of the ^-functional, i.e., 
to C(A 2 ). 
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FIG. 13: Meson masses of Fig. 1121 but in wider temperature 
region. 
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At large scales (fi > m a ) the chiral phase transition 
proceeds similar to that in the naive renormalization 
scheme jl|. In the case of the exact 0(N) symmetry, it 
proceeds through a sequence of two second-order phase 
transitions rather than a single one. In the first transi- 
tion the mean field vanishes but the meson masses still re- 
main different. The temperature of this phase transition, 
T2 — 180 MeV, is approximately /^-independent and has 
approximately the same value as that in the naive renor- 
malization scheme 0. In the second transition also the 
masses become equal, and the O(N) symmetry is com- 
pletely restored. The corresponding temperature, T cross , 
turns out to be essentially //-dependent. When the O(N) 
symmetry is explicitly violated, the sequence of two phase 
transitions is transformed into a smooth cross-over tran- 
sition. Moreover, at /1 ~ m a the results are even qualita- 
tively close to those obtained with the naive renormaliza- 
tion [l|. At small scales (say /1 < m a /2), the phase struc- 
ture becomes very complicated, however still respecting 
the NG theorem in the phase of spontaneously broken 
symmetry. 

Another result concerns both the conventional renor- 
malized HF and gHF approximations, which in fact are 
identical in the phase of restored O(N) symmetry. There 
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exists an end point of the solution, i.e., a temperature 
T enc j above which there are no solutions to the gap equa- 
tions. The occurrence of such an end point in the HF 
approximation was first pointed out by Baym and Grin- 
stein [j| and is caused by the dominant role of ln(M 2 /(J 2 ) 
terms in the gap equations, which originate from the 
renormalization procedure. The dominant behavior of 
the ln(M 2 1 a I fi 2 ) terms signals a breakdown of the HF 
approximation and the need to include higher-order cor- 
rections into the $ functional @- Another interesting 
question is whether one can find a renormalization-group 
improved <I>-derivable approximation to cure this prob- 
lem. We have found that at low scales (like /i < m a /2), 
the end point appears at rather low temperatures, leaving 
almost no room for the HF (as well as gHF) approxima- 
tion in the phase of restored O(N) symmetry. At the 
same time, at (/i > m a ), the end point moves to rather 
high energies T en( j > 1.2 GeV, hence allowing this ap- 
proximation at least at T < 1 GeV. 

Summarizing, we have found that the gHF approxima- 
tion is certainly advantageous as compared to the con- 
ventional HF one, since it respects the NG theorem and 
allows scale-independent renormalization in the vacuum. 
These properties are closely interrelated. They both re- 
quire that the set of equations of motion are degenerate in 
the phase of spontaneously broken symmetry. Neverthe- 
less, there still are serious problems with the renormal- 
ization of ^-derivable approximations for a theory with a 
spontaneously broken symmetry. At finite temperatures 
the predictions of renormalized gHF approximations es- 



sentially depend on the renormalization scale, contrary 
to the case of the renormalized perturbation theory. In 
this respect, the gHF approximation becomes similar to 
the leading-order 1/ Af-approximation, where also both 
the NG theorem and scale-independent renormalization 
in the vacuum hold true [Tfll ]. However, in view of the 
medium-independent renormalization performed in ac- 
cordance with Refs. 0, Q, 0], this scale-dependence 
at finite temperatures cannot be already interpreted as 
an artifact of temperature-dependent counter terms. It 
may turn out that this scale dependence is a consequence 
of triviality of the A</> 4 theory |25j which, when it is renor- 
malized, therefore requires an external scale to serve for 
a scale, below which it can be used as an effective field 
theory to describe the low-energy phenomenology. 
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